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Abstract: Swimming bacteria detect chemical gradients by performing temporal com- 
parisons of recent measurements of chemical concentration. These comparisons are de- 
^ ' scribed quantitatively by the chemotactic response function, which we expect to optimize 

■ chemotactic behavioral performance. We identify two independent chemotactic perfor- 

ce ' mance criteria: in the short run, a favorable response function should move bacteria up 

chemoattractant gradients, while in the long run, bacteria should aggregate at peaks of 



O . 

. chemoattractant concentration. Surprisingly, these two criteria conflict, so that when one 

' performance criterion is most favorable, the other is unfavorable. Since both types of be- 



havior are biologically relevant, we include both behaviors in a composite optimization 
that yields a response function that closely resembles experimental measurements. Our 
work suggests that the bacterial chemotactic response function can be derived from sim- 
^ ■ pie behavioral considerations, and sheds light on how the response function contributes to 

. chemotactic performance. 
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1. Introduction 

The bacterium E. coli moves up gradients to regions of high chemoattractant concentration 
by performing a biased random walk. The random walk consists of alternating runs (periods 
of forward movement) and tumbles (sudden reorientations) that arise from changes in 
flagellar rotation Q]. When the flagella rotate counter clockwise, they form a bundle 
and the bacterium swims more or less in a straight line at a roughly uniform speed. When 
one or more flagella rotate clockwise, they leave the bundle and the bacterium tumbles, 
randomly re-orienting itself |^, ^. Bacteria bias the random walk by modulating the run 
duration in response to measurements of chemoattractant concentration that are made at 
the cell surface They do not perform spatial comparisons between points along the 

cell body because of the fast diffusion across such short distances 0. 

The chemotactic response function describes how bacteria process concentration mea- 
surements to produce their behavioral run-biasing decisions. It has been measured exper- 
imentally by monitoring the rotation of single flagella on bacteria stimulated by instanta- 
neous chemoattractant pulses Q. The empirical response function is biphasic: the pulse 
provokes an immediate brief elevation of the counter clockwise probability immediately 
followed by a longer depression. We expect that the shape of the chemotactic response 
function should deliver optimal behavioral performance. 

We consider the chemotactic behavior of a bacterium at some specific position on a 
gradient of attractant. As it wanders up and down the gradient, the distribution of its 
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positions approaches a steady state. We choose performance criteria that quantitatively 
characterize the performance of the bacterium at early times in the non-steady state regime 
and at late times in steady state. Both of these regimes are biologically relevant. If the 
system navigated by the bacterium is small compared to the distance the bacterium could 
explore in the time between cell divisions (an example is bacterial aggregation into clusters 
Q), then it is the steady state behavior that matters most to the bacterium. If, however, 
the system is large - more than a few millimeters in size - or varies in time, the bacterium 
will not come to a steady state before dividing, and a single cell might never reach a 
steady state. Bacteria have no a priori knowledge of the size of their system, so their 
chemotactic strategy should benefit them in either the steady state or non-steady state 
regime. Following foraging theory |10, 11], we will assume that the chemotactic strategy 
maximizes the attractant seen by the bacterium on the timescale of bacterial divisions. 



Our first performance criterion refiects the expected velocity of bacteria at early times, 
before they have reached the boundaries of the system. It is quantified by T, a measure 
of the early time transient velocity of bacteria with a given response function. This was 
previously calculated by de Gennes [0]. Optimizing T leads to a single-lobed response 
function, which causes a population of bacteria to have a transient average velocity up 
gradients at early times. Contrary to intuition, this optimization leads to an unfavorable 
steady state distribution with bacteria accumulated in regions of low attractant. 

The second performance criterion, S, quantifies how strongly the bacteria aggregate 
about chemoattractant maxima when in steady state. Optimizing S leads to a bacterium 
that has a mean velocity down gradients at early times, but whose position distribution 
peaks at high concentrations at long times. The two performance criteria confiict: when 
one is maximal, the other is unfavorable. If both performance criteria are used to calculate 
the response function, the theoretical function closely matches the empirical biphasic bias 
curve measured by Segall, Block, and Berg |Q. The optimization procedure explains the 
curve's structure. 



Our work contributes to a body of theoretical investigations of bacterial chemotaxis. 
Schnitzer et al. [13| used Monte Carlo simulations to confirm the favorable performance 
of a biphasic response function compared to a monophasic one. Our approach supports 
their end result, though we show that aggregation can occur without a positive lobe on 
the response function. Schnitzer [14| also adopted a kinetic approach and derived results 
about steady state behavior in a variety of cases. He distinguished between 'nonadaptive 
pseudochemotaxis' and 'true adaptive chemotaxis'. In contrast, our approach emphasizes 
both transient and steady state behavior in evaluating chemotaxis. In a different approach, 
Strong et al. |15| adopted a deterministic model for tumbling and examined optimality 
in the presence of signal noise. Work by de Gennes |jl^ focused on the mean bacterial 
velocity due to a given response function. We show that this mean velocity only applies 
at early times, and we extend the framework used by de Gennes to examine steady state 
performance and performance optimization. 
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2. Model Details 



We adopt the stochastic framework used by de Gennes |12|. In this model, bacteria contin- 
uously modulate their instantaneous probability of tumbling as a function of a differential 
weighting of past measurements of chemoattractant concentration. The differential weight- 
ing constitutes the chemotactic response function, R{t). 

We assume that the chemical landscape is static and that chemoattractant concentra- 
tion is defined at every point by a function c{x). Bacteria swim along individual paths 
x{t) at a uniform speed v. The probability, P, that a bacterium tumbles in an interval 
between t and t + dt is dictated by its entire previous path, the chemical landscape, and 
the chemotactic response function: 

df \ l'^ 

P[x{t'); t]dt = — 1- dt"R{t - t")c{x{t")) (2.1) 

T' V J-oo 

where r is the mean run duration in the absence of a perturbation. In a uniform concen- 
tration, this model describes tumbling as an unbiased Poisson process with a constant rate 
of tumbling 1/f given by (1 — c J R{t)dt)/T. In a concentration gradient, P depends on 
the bacterium's history. By choosing particular forms of the response function, bacteria 
can bias their random walk so that they climb gradients and remain in regions of high c. 

We will consider first order perturbations of the Poisson process by defining R{t) oc 
a/r, with a small such that dt'R{t — t')c{t') <C 1. The constant a has units of inverse 
concentration. We expand equations as power series of such integrals and discard higher 
order terms involving products of such integrals. Equation ( |2.lD can itself be regarded 
as the first order expansion of some monotonic function of f^^dt" R{t — t")c{x{t")) that 
remains positive for all concentrations. 

In our analysis, we neglect the effects of noise due to fiuctuations of c(x). Noise 
averages to zero in all our first order expansions. The first noise contribution that does not 
average to zero is proportional to the variance of the concentration, and is of order a^c/V , 
where V is the cell volume. To neglect this term with respect to the first order term, we 
require that a <^ V . The experimental conditions described by Segall |^] correspond to 
the regime in which bacterial responses are linear and the bacteria can detect c without 
being overwhelmed by noise. 

Berg and Purcell argued that measurement integration times of about 1 second 
account for observed sensitivity to concentrations and gradients in the presence of noise. 
The response functions resulting from our analysis vary on the timescale r, about 1 second, 
so they will display biological sensitivity without explicitly requiring long integration times. 

We assume that the length scale of variations in the concentration gradient is longer 
than the average run length, so that over one run the gradient appears linear. We con- 
sider bacteria in one dimension and assume that tumbles are instantaneous and randomize 
orientation. Of course, real bacteria navigate in three dimensions and their run directions 
are not completely decorrelated by tumbles |Q, §]. Further, in real bacteria, runs directed 
up attractant gradients lengthen, while those directed downwards are the same length as 
runs in constant concentrations [Q]. Nonetheless, this simplified model gives insight into 
real bacterial behavior. 
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3. Transient Chemotaxis 



The strategic goal of a bacterium navigating a chemoattractant landscape is arguably as 
simple as producing an average velocity up the attractant gradient. De Gennes showed 
that a mean velocity can be produced when, after a tumble, a run up the gradient lasts 
longer than a run down the gradient jl^. For a population of bacteria starting at the same 
position, the expected velocity at early times will be: 

- ^ 

V ~ V O.l 

2r ^ ^ 

where At is the difference in run times moving up and down the gradient and the bars are 
averages over possible trajectories. 

The model presented in Equation ( |2.1| ) dictates that the probability of next tumbling 
at time tf after having previously tumbled at time to is 

P(t/|to) =P[x(t");i/]exp|-^*'p[x(t");t']dt'| (3.2) 

where P[x{t"); tf] is the probability of tumbling at time tf, given a path x{t"). 

Following De Gennes [12|, we consider the expected time until the next tumble: 

poo 

tixo) = / dtf{tf-to)P{tf\to) (3.3) 

J to 

where the average is taken over possible future trajectories and the bacterium is at xq at 
time to- We define t (xq) as the mean time until the next tumble for bacteria moving 
up (+) and down (— ) the gradient. We expand c(x(t)) into c{xq) ± f (Vc)(t — to), where 
V is the constant speed of a run. After expanding in a, and then using the identity 
R{t) = /q°° R{s)S{s — t)dt where 5{s — t) is a Dirac delta function, we find de Gennes's 
result that: 

V = — (t^(xo) - r(xo)) = wVVc(xo) / e-^I^R{t)dt (3.4) 
2t" Jo 

Figure |^ illustrates this integral over future paths. For bacteria to behave most favorably 
at early times, v should be maximal. 

To gain intuition about this mean velocity, consider 1000 bacteria all taking exactly 
average steps, beginning at a point xq on an infinitely long gradient, as illustrated in Figure 
H Initially, 500 move up the gradient until the time t , while 500 move down until time 
t . The average position of the bacteria is simply xq until time t , when the 500 moving 
down split into 250 moving up and 250 moving down. From this time until the up-moving 
bacteria tumble at time t"*", the mean position of the bacteria moves up at v/2. This 
phenomenon is repeated after every tumble, creating the mean velocity up the gradient 
(Equation (3J.)). In Figure |2|, this mean velocity is reflected by the thick tail of bacteria on 



the up-moving branch and the thinner tail on the down-moving branch. When the upward 
moving tail encounters a boundary on the system, bacteria are forced to tumble and the 
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Figure 1: The top figure illustrates the integration in the expression for T (Equation (3_^)). Two 
bacteria that have both just tumbled are considered as they move in different directions along the 
gradient until they tumble again at position x(tf). The bottom figure illustrates the integration in 



the expression for S (Equation (4.4)). In this case, two bacteria meet that last tumbled at points 



a;(to). One finds the expectation of their respective tumbling probabilities, P , by averaging over 
possible histories. 



mean velocity up the gradient dies away as the bacteria move toward their steady state 
distribution. Figure shows the results of a simulation that demonstrates this transient 
behavior. 



We divide out the constants in (3.4) and introduce the dimensionless performance 
measure 



-t 1 



oo 



to quantify the transient chemotactic behavior at early times. This quantity is an overlap 
integral of R{t) against a performance kernel Kq-it) = ^e"*/'^. The performance kernel 
shows the effect of the response function on the mean velocity at early times. The form 
of this kernel can be understood qualitatively. The mean velocity is proportional to the 
difference in run times between two bacteria with the same starting point, moving in 
different directions (see Figure Q) . As up- and down-moving bacteria move away from each 
other, the difference in the concentrations they measure grows until the bacteria tumble. 
Therefore a response that weights c{t) heavily in the immediate past will contribute more 
to increasing T than a weighting further in the past where concentration differences were 
smaller. This is why the performance kernel prefers recent weighting. The shape of the 
performance kernel matches simulations of the model system (Figure |3|c). The exponential 
decrease in influence of R{t) on T is due to the exponential run length of the unperturbed 
Poisson process. Note that this heuristic argument is not strongly dependent on the form 
of P chosen in Equation (2J). Any positive decreasing function of dt" R(t — t")c{x{t")) 
would yield a kernel that weights the most recent measurements most heavily. 

We can maximize T over a constrained set of response functions. We assume the 
response function to be finite and to decay to at large t. The simplest way to include 
both restrictions is to hold the integral Jq°° B?{t)dt constant. This amounts to maximizing 
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Figure 2: Solid lines indicate possible paths taken by bacteria that all execute exactly average 
paths; line thickness gives a sense of the probability weighting of each path segment. The chemoat- 
tractant gradient in this case is positive and t > t . The dotted line shows the average position 
over time: it moves to the right, indicating an expected velocity up the gradient. Note that after 
the time elapsed in this figure, more bacteria on average will have reached the furthest right point 
than the furthest left point, just because they have tumbled less frequently. 



over a set of response functions that have the same root-mean-squared deviation from 0. 
We impose the constraint 

poo 

/ R^{t)dt = o?It (3.6) 

and maximize T subject to this constraint by using a Lagrange multiplier and taking a 
functional derivative: 



T + x(^-^ Rit'fdt' 



(3.7) 



6R{t) 

Solving this condition, we calculate the optimized response function 

a 

R^i^t) = -Nrexpi-t/r} (3.8) 
r 

where Nj- is a normalization constant. This response function is proportional to the per- 
formance kernel KT{t) used to determine T; it is positive everywhere but weighted towards 
most recent times (shown in Figure^). Using this response function, bacteria moving up 
and down the gradient are progressively less and more likely to tumble, respectively. Given 
a particular tumbling position x, this results in maximally longer runs up the gradient than 



down it. A similar effect has been termed 'pseudochemotaxis' |lf:]. We call it 'transient 
chemotaxis' because unlike in pseudochemotaxis, P[x{t');t] in transient chemotaxis has a 
history dependence, and moreover we argue that short-term performance is relevant for 
bacteria in large chemical gradients. 

Surprisingly, although Rr{t) maximizes the expected velocity up the gradient, it leads 
to an unfavorable steady state distribution in which the bacteria spend more time in low 
chemoattractant concentration regions. The simulation in Figure |3|b shows the initial 
favorable transient velocity and the unfavorable steady state for an all-positive response 
function. This counter-intuitive result can be explained as follows. Imagine two bacteria 
passing each other on a linear concentration gradient (see Figure |l|) . The one heading down 
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the gradient has high c in its past, so its value of / Rr{t — t')c{x{t'))dt' is larger on average 
than that of an upward-moving one at the same position. Equation ( |2.1[ ) then shows that 
the probability of tumbling is lower for the bacterium moving down the gradient. Since this 
is true at all points on the gradient, more bacteria will accumulate in the low concentration 
areas. The unfavorable steady state of a positive response function was previously shown 
in numerical simulations [^] and noted in Schnitzer's analysis p^ . 

Initial velocity need not indicate the eventual steady state distribution, as the following 



thought experiment shows (discussed by Lapidus Schnitzer et al. [13|, and Schnitzer 
fl^ ). Consider a closed tube containing a gradient in density of steel wool. At one end of 
the tube, mean free paths of a molecule are short, while at the other end they are long. 
After each collision, because of the gradient in the wool, a molecule has an expected net 
displacement towards the sparse end of the tube. In steady state, however, gas molecules 
are distributed evenly throughout the free volume of the tube. Therefore, although the 
expected net displacement after each collision creates an initial mean velocity towards 
the sparse end, it does not determine the steady state distribution. For gas molecules, the 
collision probability is determined by a particle's instantaneous position. For bacteria using 
an all-positive response function, both and t are longer in higher c regions because 
Jq°° Rr{t)dt / 0, making path length depend on position. It is the history-dependence of 
Rr that causes the bacteria to aggregate in regions of low c. 



4. Steady State Bacterial Distribution 

Here we show how the steady state distribution of bacteria depends on expected tumbling 
rates. The expected tumbling rate for a bacterium at position x depends on whether it is 
moving up or down the gradient, and is given by P (x) or P (x), respectively, where bars 
are averages over possible histories ending at x. In steady state these averages will not be 
functions of t. 

In steady state, bacterial flux is zero and the bacterial steady state concentration, h{x), 
can be written in terms of the probabilities P (x) (see the Supplementary Information): 

5(x) = 60 exp |£ ^ (p-(x') - P+(x')) } (4.1) 

This reproduces a more general result derived in [^]. 

With net flux zero, the number of upward-moving bacteria must equal the number of 
downward-moving bacteria at any point x. If P (x) 7^ P (x), then the fraction of bacteria 
passing through a point from the left will not equal that passing through from the right. 
In order to keep the number fluxes equal, the number of bacteria on each side of that 
point must be different. Maintaining this balance generates the form of the distribution 
in Equation ( [4.1| ). When the tumbling rate is higher for down-moving bacteria arriving at 
point X, bacteria aggregate at the top of the gradient in steady state. 
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Figure 3: We performed discrete time simulations of the model on a positive concentration gradient 
with reflective boundary conditions to see the result of different R{t) on transient and steady state 
behaviors. Bacteria were released from the center of the gradient (a, left) and evolved until they 
arrived at a steady state distribution. R{t) was chosen to weight positively only at 6 seconds before 
the current time, t (a, right); that is, it weights only c(t — 6). It was further chosen so that the 
maximum perturbation from the average tumbling probability was 30%. b In a gradient of length 
60uf, bacterial distributions and the mean position of bacteria were found using a response function 
with 6 = f, where f is the run duration averaged over the box. At early times, bacteria are clustered 
and have a mean velocity up the gradient. After the bacteria hit the boundary, they approach a 
steady state peaked at low c. Note that more bacteria have reached the righthand wall than left 
hand wall at t = lOOf . For this response function, T > and <S < 0; both results are reflected in 
the bacterial behavior, c We varied 9 and calculated T from the initial slope of the lower plot in b. 
The result shows the contribution of R{9) to T. The solid line is the transient performance kernel, 
Kr, derived in the text, d In a short length scale gradient (4uf), we varied 6 and calculated <S 
from the bacterial distributions at long times. This shows the contribution of R{0) to <S. The solid 
line is the steady state performance kernel, Ks, derived for a similar situation (see Supplementary 
Information). Error bars in c and d are 1 SEM. 
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We now express the tumbling probabilities in terms of R{t). In order to calculate 
P (x) we must consider all possible histories of bacteria reaching point x. Histories and 
instantaneous tumbling probabilities both depend on R{t), so the difference P (x) — P (x) 
that governs steady state aggregation can be expressed in terms of the response function. 
By integrating over paths for bacteria arriving at x (details of the derivation are in the 
Supplementary Information), we find that 

P"(x) - P^{x) = 2vVc{x) / -{t/T + t^/2T'^)e-*/^R{t)dt (4.2) 

Jo 

This integral should be positive to obtain an advantageous steady state distribution 
with more bacteria at high concentrations. The x dependence in Equation ( |4.2| ) comes 
through the Vc factor, which integrates immediately to c(x), giving the steady state dis- 
tribution: 

.(.)=*„»p{.(.)(^:|^^)} (4.3) 

The quantity in round brackets does not depend on x. We introduce the dimensionless 
version of this quantity, 

= ^^WT = I [ -"/^ ^ 'V2r>-'/^H(,)d, (4,4) 

as a performance measure of the steady state distribution. This is an overlap integral, with 
a performance kernel Ks{t) = ~^{t/T + t"^ /2T'^)e~^^'^ . A response with large S yields a 
steady state distribution with the bacteria aggregated favorably in high c regions. 

When S is maximized by the same procedure used in Equation (|3.7| ), one finds a 
response function 

Rs{t) = --Ns{t/T + tV2r') exp{-t/T} (4.5) 
r 

which is negative everywhere, zero at t = and at long times, and peaked at t = 
(see Figure ^). The negative values of this response function mean that bacteria moving 
down the gradient at point x, with high concentrations in their past, have higher tumbling 
probabilities than bacteria moving up the gradient at x, with lower concentrations in their 
past. 

Because Rs{t) is negative, it results in t {x) > t^{x), creating a transient velocity 
down the gradient at early times. Although this response function gives a beneficial steady 
state distribution, it yields detrimental behavior at early times. 

One can understand the steady state performance kernel qualitatively. The perfor- 
mance measure S considers the difference in tumbling probability between two bacteria 
at the same point in space but coming from opposite directions (see Figure In this 
case, measurements of c are most different in the past, while the most recent concentration 
measurement, c(x), is the same for both bacteria. This weighting is reflected in the perfor- 
mance kernel Kg and in the optimal response Rs{t), in which concentrations in the past 
are more heavily weighted. Concentration measurements in the more distant past could 
have been made where Vc was different from the current Vc and cannot be reliably re- 
lated to the current gradient. Therefore, such distant information is not useful for making 
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run-biasing decisions and is not weighted heavily by the kernel |17]. Figure ^ shows the 
derived performance kernel and results of simulations of the model in a small system. 
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Figure 4: a Response functions that optimize the performance measures T, S, and T + AS, where 
A = 1/2. Note that all three functions are normalized such that J R{t)^dt = o? It. b The points are 
data from Figure 1 in |^ showing the counter clockwise bias in flagellar motor rotation after a very 
short impulse of chemoattractant at time i = 0. The bias response is linear in this experiment's 
regime. The solid line is a best fit of Rs,T{t) to the data, using a 10 Hz low-pass Gaussian filter 
to realistically smooth discontinuities. The fitting parameters were A, r, and an overall amplitude, 
and the least squares fit was A = 0.56 and r ~ 0.9 seconds. The bias of a single fiagellum is related 
to the tumbling probability P[x{t');t\, but is not identical because multiple flagella are involved in 
running and tumbling Q and cooperative effects could be involved. 

5. Optimizing the Response Function 

The response functions resulting from optimizing the two performance criteria have oppo- 
site signs, so that optimizing T leads to an unfavorable S and vice versa. Both aspects 
of performance are biologically relevant - bacteria should move up gradients when not in 
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steady state and remain at high concentrations as they approach steady state. We expect 
bacteria to optimize a composite criterion that preserves both aspects of performance. One 
can imagine a variety of ways to maximize a combination of the two quantities, but maxi- 
mizing any positive increasing function of both T and S will produce a solution that is a 
linear combination of Rg and Rt- We therefore adopt the most straight forward way and 
maximize the quantity 

T[R{t)] + A S[R{t)] (5.1) 

where A is some unknown weighting factor of the two performance measures. As before, 
we constrain R{t) and take the functional derivative of this equation to find a response 
function that compromises between maximizing the expected run length up gradients and 
the steady state bacterial distribution. That response function is 

RrAt) = ^Nr,s exp{-t/r} (1 - A{t/T + tV2r')) (5-2) 

which is proportional to Rr + ARs . 

It is reasonable to set J R = because there are physical bounds placed on the run 
length of real bacteria. Purcell |18| pointed out that run duration should be chosen at least 
large enough so that, for a given v, a bacterium outruns the diffusion of the chemoattractant 
c during its run. This lower bound on run duration is roughly 1 second. Further, in real 
situations runs longer than roughly 10s are turned 90 degrees off-course by rotational 
diffusion [Q, |l^, setting a maximum useful run duration. Neither of these limits depends 
on c. Bacteria should be sensitive to gradients by maintaining a large a but must not 
allow their run durations to wander outside these bounds in widely varying concentrations. 
Run duration is dependent on a Vc term as derived in the text and on c J" R{t)dt. The 
integral of R should be zero to allow for sensitivity to Vc while keeping f within the limits 
above, thus creating a large dynamic range for the response ||l^. This argument leads us 
to set A = 1/2 so that J R = 0. Experimentally, Alon et al. [20| have shown that this 
robustness of run duration to changes in absolute concentration is a property of the E. coli 
chemotactic network when cells respond to aspartate (though see [||). 

The optimized response function is shown in Figure^, where we have required f R = 
0. It predicts a sharp positive immediate response with a drawn out negative response 
peaking around 2.5r. It was obtained here by developing a response function that (1) 
maximizes the transient velocity of bacteria up gradients when they are not in steady state 
and that (2) creates a steady state where bacteria aggregate in high concentration regions. 
The initial, short, positive lobe in RT,si't) makes T > and serves to move the bacteria 
up gradients when not in steady state, while the second, longer, negative lobe makes 5 > 
and serves to produce the advantageous steady state distribution. 

The functional form of Equation ( |5.2| ) fits the actual response function exhibited by 
individual flagellar motors in Segall et al. |^] (see Figure ^b). Our theory concerns the 
whole cell, not single flagellar motors; correlations between the activity of single flagella 
and the behavior of the whole bacterium are not well understood [21, Nevertheless, we 
find a surprisingly good fit. We have left A and r as fitting parameters, and the best fit 
yields A = 0.56, which matches our expectation that A ~ 0.5. 



- 11 - 



Values of T and S for any response function can be easily found by calculating their 
overlap with the kernels in Equations { \i.5\) and ( [4. 4]) ; a summary of such calculations 
for our three optimizations is shown in Table 1. The top half of the table provides the 
qualitative picture independent of model details, while the lower half provides the values 
of T and S given by our model. 

Table 1: Transient velocity values (T) and steady state strength of aggregation (S) for the 
various response functions R(t). The first set is heuristic, derived from qualitative arguments, 
while the second set is derived from our particular model. The first item in each section is the 
response function maximizing T, while the second maximizes S. The third maximizes both T and 
S, as described in the text. More positive values of T and S indicate more favorable behavioral 
performance; in the heuristic section, favorable values are represented by (+) and unfavorable by 
(-)■ 



Response Function 


Equation 


T 


S 




Positive lobe, weighted towards t = 


+ 




Rs{t) 


Negative lobe, weighted towards t ~ r 




+ 




Initial brief positive lobe; negative lobe peaked beyond r 


+ 


+ 


Rr{t) 




0.7 


-0.5 


Rs{t) 


-^Ns (t/r + tV2r2) e^'^^ 


-0.4 


0.9 


RT,s{t) 


^Nr,s{l-\{t/T + t^/2T^))e-''l- 


0.5 


0.05 



6. Discussion 

The biphasic shape of the chemotactic response function has long been interpreted as a 
temporal comparator that automatically adapts to offsets in chemical concentration 
Here we have examined two aspects of bacterial chemotactic performance and found that 
neither aspect optimized alone produces a biphasic response. A composite response func- 
tion, simultaneously optimizing both measures of performance, closely fits the shape of the 
experimental data. This leads to an additional interpretation of R{t) as optimized with 
respect to these two behaviors, connecting each lobe to distinct behavioral performance. 

Our theory makes novel predictions about the behavior of wild-type and mutant bac- 
teria. The functional fit of i?T,5 to the wild-type data is quite good, so that we predict that 
experimental measures of wild-type T and S in the linear regime should roughly match 
those in Table 1. If L is the decay length of b[x) when bacteria are on a linear gradi- 
ent, it is related to the expected transient velocity up the gradient, Vd, by the relation 
L = [v'^r/vd) * (T/S). The first factor could be found on dimensional grounds, but we 
predict T/S ~ 14 for wild-type E. coli. For a gradient that elicits Vd = ^^J^-Js, this predicts 
L ~ 5 millimeters. Response functions of mutant bacteria can be calculated (see [^2| or 
||2^ ) or measured experimentally, as Segall et al. |^] have done for strains with mutant 
cheZ and for strains with mutant cheRcheB. Both these mutant response functions are 
entirely positive with durations of roughly 5 and 1 seconds, respectively. We predict that 
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both mutants will have transient velocities up gradients, but that both will reach an unfa- 
vorable steady state distribution. Available data for both mutants does not rule out these 
predictions [24, 25, 26]. Microscopic observations oft^ — t or measurements of 6(x) in 
static spatial gradients could evaluate the validity and limits of this theory. 
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A. Steady State and Tumbling Probabilities 

Here we derive Equation ([4.1|), the steady state distribution of bacteria in terms of the 
expected tumbling probabilities of bacteria coming from the right and left. 

In the steady state, the net flux of bacteria through a point between x and dx must be 
zero. Bacterial concentration at point x is represented by b^{x), where the (it) indicates 
the direction of movement. Setting the net flux to zero yields 

b+{x-dx){l- -P^ix)—) + b-{x-dx)-P (x)— = (A.l) 
2 V 2 V 

b-{x + dx){l - -p-{x)—) +b+{x + dx)-P'^{x) — 
2 V 2 V 

The left hand side counts bacteria moving up the gradient past x, while the right hand side 

counts bacteria moving down the gradient past x. The probability of actually reversing 

directions is half the probability of tumbling. The first term on each side represents bacteria 

continuing on their present course; the second represents bacteria passing through the point 

after an instantaneous reorientation. Retaining only zeroth order terms gives b'^{x) = 

b~{x), so one can replace each by b(x)/2. Expanding b{x) about x and retaining only terms 

up to first order in dx yields a differential equation governing the steady state distribution 

of bacteria: 

Vbjx) _ p-{x)-P^{x) 

This is integrated to find Equation ( [4.1D . The zero flux condition holds at all points in the 
system. The equation is modified at system boundaries, depending on conditions there, 



but the first order differential Equation ( A. 2 ) describes the bacterial distribution far from 
the boundaries. 



B. Optimizing the Steady State Distribution 

Here we calculate P {x) — P^{x) in the steady state and find a response function that 
optimizes the bacterial distribution, b{x). 
We are interested in: 

P~(x) -P+(x) = - / dt'R{-t')[c+{t')-c-it')] (B.l) 

T J-oo 
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where superscripts indicate that the averages are taken over all paths ending at the position 
X moving either up (+) or down (— ) the gradient at x. We take the gradient to be positive, 
so that up-moving bacteria are right-moving bacteria. 

To calculate this quantity, we must insert averages over possible paths. Assuming that 
the most recent tumble occurred at time to) the one before that at time ii, and so on, we 
will average over the values of toi ^ii • • • and over the directions of the runs during each of 
the intervals. Essentially, because is multiplied by R{t) in the integral, this calculation 
can neglect first order effects of R{t) on c^: runs become simply exponential in length and 
there is an equal probability that a bacterium came from left or right before its last tumble. 

We first average over the directions of the runs. In the steady state, for a tumble at tj, 
the probability Q~^{ti) that the bacterium was moving up the gradient before the tumble 
will be proportional to the population of bacteria found to the left of x{ti). That is 

Q+(^ N ^ b{x{ti) - vjU - tj+l)) 

b{x(ti) - V{ti - ti+i)) + b{x{t^ + v{U - U+i)) ^ 

Now assume that b{x) varies slowly over one run and expand b{x) to first order in 
ti - U+i. Then 

b{x{ti))-VbviU-ti+i) If Vb \ 

Inserting just the averages over the directions of motion in the various intervals and 
leaving the ti fixed for the moment, we can expand Equation (|B.l| ) as: 

j\dt'R{-trc^{t')^ (B.4) 
!ldt'R{-t')cHt') 
+ //; dt'R{-t') [Q+ (to)c(±+) {t') + Q- (to)c(±-) (tO] 
+ //;dfi?(-t')[Q+(to)Q+(ti)c(±++)(t') + Q+(io)Q-(ti)c(±+-HO 

+g-(to)Q+(ti)c(±-+Hi') + Q'{to)Q-{tM^--\t')] + ... 

Here c^"' (t), for instance, denotes the concentration seen by a bacterium that has moved 
in the -|- direction after to, moved in the -|- direction in the interval [ti,to]) and in the 
— direction in the interval [t2 5ii]- From Equation ( |A.2| ), Vb/b = {P — P )/v, which 
is a sum of terms proportional to integrals of the form J R{t — t')c{t')dt' , so we may set 
Q^{ti) = 1/2 to keep only terms to first order in f R{t — t')c{t')dt' . 

Now we look at the averages in square brackets in Equation (B.4). In the first such 
average, we expand c{t) about to find that 

Q+(to)c(±+)(t') + Q-(to)c(±-)(i') (B.5) 
= l/2(c±(to) + v{t' - to)Vc + c±(to) - v{t' - to)Vc) 
= c±(to) 



The second square bracket in (B^)reduces to {c^^~^\ti) + c*-^ \ti))/2 and the rest of the 



square brackets will reduce similarly. 
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Inserting the averages over the tumbhng times tQ,ti, . . ., Equation (BJ) expands to: 

f_^dt'R{-trcHt')^ (B.6) 
!l^dtoDito\0)f,yt'Ri-t')c^it') 
+ dtoD{to\0) ft'^ dhD{h\to) f!^ dt'R{-t')c^{to) 
+ dtoD{to\0) ft'^ dhD{h\to) dt2D{t2\h) 

where D{62\6i) is the probabihty that the bacterium tumbled at 62 given that it tumbled 
later at 9i, where 62 < 9i so that we are reconstructing the tumbles backwards in time. 

It remains to write out the factors D{ti^i\ti) explicitly in terms of R{t). Given the 
model, 



f U 
exp<^ - / dt'P[x{t");t' 

/ dti+iexpj- / dt'P[x{t");t' 
-00 I ti+i 



where the expression is normalized to integrate to 1. Keeping only terms up to first order 
in R{t) allows us to set D{ti+i\ti) = 1/r exp{(ti_|_i — tj)/r}. Then P (x) — P (x) becomes 

P~{x)-P^{x)= (B.8) 
/-oo ^e*o/^ /° f i?(-t')[c+(tO - c-{t')] 

+ /-oo /-oo T^e*^/^ /? '^R{-t')[c-^{t0) - C-{to)] 

1 fO dto r*o dti rh dt2 to. It 
' J— 00 r J— 00 T J — 00 r 

X M + + )(fl) + c( + -)(fl)-c(- + )(tl)-c( — ^ 

First consider the quantities in square brackets. The first one can be approximated 
by c^{t') — c~{t') = 2vt'Vc. The second one can also be approximated the same way, but 
tQ replaces t' . Whether the third one can be made proportional to the gradient of the 
concentration depends on how quickly the gradient varies in space. Several of the terms in 
(B.8) might be approximated in terms of Vc; it is a question of what sorts of gradients a 



bacterium typically encounters. 

Consider a generic chemical landscape that is approximately flat on a large length 
scale L. On length scales smaller than L the concentration may vary significantly. In such 
a landscape, quantities like those in square brackets in (|B.8| ), representing measurements 
made in the distant past ^ will have a tendency to sum to zero. In particular, 
if L ~ 2vT, the quantity in the third square bracket in (|B.8| ) will be close to zero since on 
average 



c(++)(ti)+c(+-)(ti) ~c(-+)(ti)+c(— (B.9) 

Terms coming from the past where \t'\ ^> cease to contribute to P [x) — P^ {x). 
Where we cut off the series in Equation ( [B.8[ ) is a biological question. We expect bacterial 
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strategy to make minimal assumptions about the extent of the gradient, and therefore we 
will cut off the series after only a few terms. We drop all terms that refer to times earlier 
than ti. The resulting expression for P (x) — P (x) has two terms proportional to Vc. 
The terms proportional to Vc are precisely the ones that allow us to optimize the steady 
state distribution in a concentration independent way. 
We find that 

P-{X)-P^{x) = ^-i^[f_^dtoe'^/^J,yt'R{-t')t' (B.IO) 
+ f-oodtot'^dt,e^^/rf^yjLR^-t')to] 

Letting R{t) = dsR{s)5{t — s), we find the simple result that P (x) — P~^{x) can be 



written as the overlap in Equation (4.2) 



Keeping more terms in Equation (|B.S| ) corresponds to assuming that the bacteria 
propagate on gradients with variations of longer length scale. Such an assumption produces 
a kernel of e"*/^ multiplied by an expansion of (1— e*/"^) in t/r. Retaining an infinite number 
of terms leads to a kernel proportional to e~*/^ — 1. As long as a finite number of terms are 
retained, the qualitative features of the performance kernel don't change: it starts at zero, 
peaks negatively, and returns to zero for large t. If the number of retained terms remains 
small (less than 5, for example), the quantitative features are not much affected either. 

Figure ^ illustrates this discussion. It shows the results of simulations of the model 
that demonstrate the contribution to S from each portion of the response function. The 
response functions are chosen as in Figure ^ to weight only c{t — 6) in determining the 
turning probability. Bacteria navigating concentration gradients of various length scales 
are considered. On the timescales shown, the white points in the figure represent bacteria 
in an effectively infinite linear gradient. Bacteria heading down such a gradient at x on 
average have a higher c history for all previous times than ones heading up the gradient 
at X. Therefore, measurements made any time in the past, including long ago, affect S 
by the same amount. On the other hand, when the gradient is not infinite, as in any real 
case, Vc will change on some length scale. Measurements of c{t) made with a delay large 
compared to the time to traverse this length scale will average out in the sum over histories. 
Such measurements therefore cease to affect S. As the grey points in the figure show, on 
a relatively short length scale, measurements of concentration made 20f in the past no 
longer contribute to 5. In the time of 20f, the bacterium is likely to have bounced against 
a wall (or moved over a peak, if one views the reflecting box as an infinite triangle wave), 
so such measurements no longer reflect the bacterium's current gradient. With still shorter 
gradients, such effects become more pronounced. The black points in the figure are on a 
length scale that roughly matches our theory, in which we restricted bacteria to looking 
at only the previous two runs. The simulation conditions mean that bacteria cannot use 
measurements made long ago, while our theory posits that bacteria should make minimal 
assumptions about gradient length. Both lead to similar kernels showing the infiuence of 
R{e) on S. 
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Figure 5: Using the same parameters as in Figure l^d, we here show the results on different length 
scale gradients. Such gradients with reflective boundary conditions can also be thought of as infinite 
triangle waves, a more natural picture than a box. Black points are simulations of bacteria in a 
box of total length 4uf , grey points are in a box of length 8wf , and white points are in a box of 
length 80wf . The dotted curve shows the performance kernel for S in the case where the bacterium 
assumes a linear gradient throughout its entire history (an infinite gradient). The solid curve is the 
performance kernel that makes minimal assumptions about the extent of the linear gradient. 
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